Shear viscosity and spectral function of the quark matter 
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We discuss the shear viscosity of the quark matter by using the Kubo-Mori formula. It is found 
that the shear viscosity is expressed in terms of the quark spectral function. If the spectral function 
is approximated by a modified Bright- Wigner type, the viscosity decreases as the width of the 
spectral function increases. We also discuss dependence of the shear viscosity on the temperature 
and the density. 
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The existence of the quark gluon plasma (QGP), which 
is predicted by Quantum chromodynamics, has not been 
discovered in Nature. In order to produce such a new 
state of matter, the experimental investigation started at 
the Relativistic Heavy Ion Collider (RHIC). The data at 
RHIC, however, seems to reveal some unexpected prop- 
erties of the high density matter produced in the exper- 
iment 0-0. Namely it could be explained by a fluid 
model with small viscosity; it is almost perfect fluid. This 
fact encourages many researchers to calculate the trans- 
port coefRcients of the quark matter 5]- 20]. Moreover 
man y st udies by the lattice QCD have been in progress 

Mm 

It is the purpose of this note to calculate the shear vis- 
cosity of the quark matter. There are two kinds of par- 
ticles in the QGP: quarks and gluons. In this paper, we 
take up only the quark sector so that the Nambu-Jona- 
Lasinio (NJL) model is available for the quark matter. 
As for calculation method of the viscosity, we use the 
correlation function method which is called Kubo-Mori 
formula ll-lll. 

According to the Kubo-Mori formula 27J, the shear 
viscosity 77 (T) at temperature T is given by 



dte'' 



dr(J,,(r,t), J,^(0,0)), (1) 



where J^y is the x, y component of the energy-momentum 
tensor of the quark matter. The correlation function in 
the right-hand side is defined by 
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(2) 



where A and B are operators of any physical quan- 
tities and H denotes our Hamiltonian. The bracket, 
{A) = TT{Ae~^^)/Tie~^^ , means the thermal average 
at temperature T {(3 = 1/T). Using partial integration in 
the right-hand side of Eq.(l), the viscosity can be trans- 
formed into 



r7(a;) = -[n»-n«(0)]. 

UJ 



(3) 
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Here 11-'^ (a;) is a retarded Green's function defined by 

n«(c.) = -z j^"dte^"*y dr([J,,(r,i),J,,(0,0)]), (4) 

where [ , ] in the integrand denotes the commutation 
relation. Noting that (n^(w))* = n^(-w), the (static) 
viscosity is reduced to 



r] = ri{uj = 0)= — — Imn^(w) 



(5) 



UJ = + 



In order to calculate the above 11^ (cj), it is convenient to 
transform into the imaginary time (Matsubara) formal- 
ism. We introduce the following correlation function, 

n(a;„) ^ - dre-'^-^ j dr(r,(J,j,(r, r) J,j,(0, 0))), 

(6) 

where the Matsubara frequency is represented by utn — 
2nnT and Tr means the (imaginary) time ordering op- 
erator. As is well known, the retarded Green's func- 
tion n^(w) is obtained by the analytic continuation: 

n«(^) = n(a;„)|,„^_^,,. 

We take the NJL model for the quark matter in this 
paper j30l | . Then the canonical energy- momentum tensor 
is read as 



(7) 
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where ^ is the field operator for quarks. If this expression 
is substituted into Eq.(6), the correlation function 11 is 
written as 

nK) = ^E / (^P.'Tr[72G(p,c.z+c.„)7'G(p,c.O], 

(8) 

where G{p, uj) is the full propagator for the quark. This 
expression corresponds to the Feynman diagram drawn in 
Fig.l. Of course there are many higher-order diagrams 
involved in this process, such as the ring (bubble) dia- 
grams like the Nambu-Goldstone (pion) mode [sj- But 
in the present case, such terms do not contribute to the 11 
because the trace factor in the 11 is suppressed due to the 
presence of the 7^. Hereafter we are restricted ourselves 
to the simple diagram of Fig.l. 
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it is rewritten as 
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FIG. 1: The Feynman diagram for the correlation function 11 
in Eq.(8). 
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FIG. 2: The countour for the calculating the integral of z in 
Eq.(lO). 



Now let us follow the procedure taken for the calcula- 
tion of the electrical conductivity in R ef.fs^. We consider 
the spectral representation for the full propagator, which 
is written as 



r°° de 

Ga/3{p,UJi) = / — 
J-oo 



de Pa/3(P,£) 



lUJl 



(9) 



We substitute this expression into the correlation func- 
tion (Eq.(8)) and the summation over Matsubara fre- 
quency is replaced by the contour integral: 



dz 
27ri 



n{z)TY[G{p, 2)7'G(p, z + ic^„)7'](10) 



where n{z) = (1 + e^'^)^^ is the Fermi distribution func- 
tion. The countour G is divided into three pieces as 
shown in Fig. 2, because the integrand has branch cuts 
on the two lines z = s and z = e — iw„ where e is real 
(Note that poles of n{z) do not lie on the two branch 
cuts). Since the integral along the large circle vanishes, 



ds_ 
2Td 



n(e)Tr[G(e + iS)-f^G{e + iw„)7^ 



-G{e - i6)j^G{e + iujnh^ + G{e ~ iw„)7^ 
xG(£ ^ ^5)7^ - G(£ - zw„)7'G(£ - 

where b is an infinitesimal positive number introduced 
in order to avoid the branch cuts. Noting the relation, 
G(e -|- ib) — G(e — iS) = —ip{e), the above equation is 
expressed as follows: 



d£ 
27r 



n(£)Tr[(G(e + iujn) + G(e - iLJn)h^p{eh^ 



(12) 

If the analytic continuation icun — > LU + iS is carried out, 
the imaginary part of the above equation is 



Im5 = 



de 1 
27r2 



{n{e + uj)- n(£))Tr[p(£ + uj)-f^p{e)j^ 



(13) 

because Tr[p7^p7^] is real. As a result we obtain the 
shear viscosity expressed by the quark spectral function, 
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de 
2^ 



^pJ^TMeh'piem (14) 



This equation means that the calculation of the shear 
viscosity is reduced to that of the spectral function. The 
appearance of dn/de represents the Pauli blocking effect 
(See Eq.(18)). 

Next stage is to discuss the quark spectral function. 
We take the following simple parameterization for the 
retarded (advanced) Green's function. 



G«(p,c.) = 
G^{p,uj) = 



1 



p • 7 - M + isgn(po)r ' 
1 

p • 7 - Af - isgn(po)r ' 



(15) 



with the definition of po = uj + p, {p: the chemical po- 
tential). Here M and F represent the effective mass and 
the width of the quark respectively. These quantities are 
not determined and regarded as parameters in this pa- 
per. From these equations, we get a spectral function of 
modified Bright- Wigner type. 



p(p,w) = i(G^(p,^ 



G^(p,c.)) 



2Fsgn(po) 



(p • 7 - M)2 + r2 ■ 

(16) 

It should be noted that there exists a sum rule 
/ p{p,uj)duj/2TT = 7° for the spectral function. This 
equation is satisfied approximately when the value of F 
is small: F < M. If this expression is substituted into 
the trace in Eq.(14), it is written as 

Tr[p7V7'] = ^^^(8AfV, - X)T\ (17) 

Here and Nf denote the numbers of the color and 
flavor of the quark and the use is made of the definition: 
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FIG. 3: The shear viscosity as a function of the width F at the 
temperature T = 150MeV (brack circles) and T = 200MeV 
(white circles) in the case of /tt = lOMeV. 
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The shear viscosity as a 


function of width 



the chemical potential /i — lOMeV (black circles) and /i = 
lOOMeV (triangles) in the case of T = 150MeV. 

X = (p2-M2+r2)2+4M2r2. Substituting this equation 
into Eq.(14), the shear viscosity is written as 

^ T y 277 7 (27r)3 [(p2_M2+^2)2+4M2^2]2• 

(18) 

where we have neglected the second term in Eq.(17) be- 
cause the main contribution in the above integral comes 
from the integral region where X is small. 

Now let us discuss the results of the numerical calcula- 
tions of the shear viscosity. Since we are interested in the 
QGP phase, our calculations are restricted to the case of 
high temperature or high density. As for the effective 



quark mass, wc take M = lOOMcV {Nc = 3 and Nf = 2) 
for convenience. In Fig. 3, the shear viscosity is drawn as a 
function of the width F at T = f 50MeV and T = 200MeV 
in the case of = lOMeV. Similarly Fig. 4 shows the 
shear viscosity as a function of T a,t /j. — lOMeV and 
1^ = lOOMeV in the case of T = 150MeV. These figures 
evidently show that the viscosity is a rapidly decreasing 
function of the width F. The viscosity seems to diverge at 
F = where the quark matter becomes ideal gas. This is 
understood by noting that (r/(a;^ -|-r2))2 (ttS^x))^ as 
F — > 0. As the r becomes larger, the viscosity decreases 
rapidly. This is also apparent because (F/(.t^+F2))2 Q 
asT oo. The width means the inverse of the quark life 
time and reflects the strength of the interaction. There- 
fore the quark matter becomes almost perfect fluid un- 
der strongly interacting state. This property is consistent 
with the recent RHIC data. But we do not discuss the nu- 
merical values of the shear viscosity in comparison with 
the experimental data. The reason is that the spectral 
function used in this paper is a phenomenological one 
and has not been determined theoretically. 

We also comment the dependence of the shear viscosity 
on the temperature and the density. From Figs. 3 and 4, 
it is seen that the viscosity increases as the temperature 
increases. As for the density, it is a slowly increasing 
function. These dependences seem to be consistent with 
the classical expression of the viscosity of gas: rj = pvl/3 
(p=density, t;=velocity and 1= mean free path of the par- 
ticle). 

In conclusion, we have calculated the shear viscosity of 
the quark matter with the use of the Kubo-Mori formula. 
We have obtained the shear viscosity expressed in terms 
of the quark spectral function. Assuming a modified 
Bright- Wigner type of the spectral function, the shear 
viscosity is a rapidly decreasing function of the width; the 
quark matter approaches to perfect fluid under strongly 
interacting state. 
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